On the Nature of the Rapidity- Spectra at RHIC and Some Other 

Energies. 



Goutam Sau""^* S. K. Biswas^jA. C. Das Ghosh^f A. Bhattacharya'^lc S. Bhattacharyya^^ 

^ Beramara RamChandrapur High School, 
South 24-Pgs,743609(WB),India. 
^ West Kodaha Adarsha Siksha Sadan, 
New Barrackpore, Kolkata-700131, India. 
^ Department of Microbiology, 
Surendranath College, Kolkata- 700009, India. 

^ Department of Physics, 
Jadavpur University, Kolkata- 700032, India. 
^ Physics and Applied Mathematics Unit(PAMU), 
Indian Statistical Institute, Kolkata - 700108, India. 



Abstract 

On the basis of the Grand Combinational Model (GCM) outlined and somewhat detailed in 
the text, we have attempted to capture here the several interesting assorted characteristics of 
the rapidity-spectra of the major varieties of secondaries produced in diverse nuclear reactions 
at various energies, though the main thrust of our work lies on addressing the data-trends from 
RHIC-BNL experiments. Obviously the core of the present approach is purely phenomenological. 
Still, the method and the model address the features of the data modestly well. And the 
method appears to have the rich potentiality, if the systematic sets of data for rapidity-studies 
at gradually increasing energies were available. 

Keywords: Relativistic heavy ion collisions; inclusive cross-section. 

PACS nos.:25.75.-q, 13.60.Hb 



* e-mail 
^e-mail 
■'e-mail 
^e-mail 
^e-mail 



sau_goutam@yahoo.com 
suniLbiswas2004@yahoo.com 
dasghosh@yahoo.co.in 
pampa@phys.jdvu.ac.in 

bsubrata@www.isical.ac.in (Communicating Author). 



1 



1 Introduction 



Unlike the cases of pT-spectra, theoretical (even phenomenological) , studies on the rapidity-spectra 
are relatively much fewer, though these two observables [{pt^y) or (pt^v)] just twins for ana- 
lyzing the nature of invariant cross-sections for the secondary particles produced in the arranged 
laboratory collisions at very high energies. In fact, the experimentalists also are not immune from 
such biases, mainly due to the accompanying difficulties in the measurement processes. And the 
net effect of all these put together is the number of data sets on rapidity spectra is obviously much 
less abundant and they suffer from large measuremental uncertainties or large errors of both sys- 
tematic and statistical nature. Despite these strong and somewhat insurmountable limitations, in 
the present work we will confine ourselves to the studies on rapidity-spectra produced in Pb-|-Pb 
reactions at marginally high energies and in RHIC-BNL interactions involving Au-|-Au collisions 
at ^/sNN = 200 GeV in the light of Grand Combinational Model (GCM). The stimulus stems 
from the fact that in the past we analyzed a large part of the data on p^-spectra for RHIC-BNL 
collisions and a considerable part of the rapidity-spectra at CERN SPS-energies. So, for the sake 
of completeness of our studies on the behaviour of this specific model (GCM), we need to check 
whether this model as such (or any extension of the model for higher energies) could explain the 
relatively sparse data on the rapidity-spectra and some other related properties involving these 
experimental measurements on rapidity-dependent physical observables. 

The plan of this paper is as follows. In Section 2, we attempt to provide an outline of the Grand 
Combinational Method (GCM) and give the successive procedural steps. The Section 3 gives an 
inkling into the basic root of the expressions used in Section2. The Section 4 exhibits the results 
and some discussions on them. The last Section (Section 4) contains the final summary and the 
conclusive remarks. 

2 The Working Model : An Outlook 

Following Faessler[T], Peitzmann[2] and also the work of Schmidt and Schukraft[3], we propose 
here a generalized empirical relationship between the inclusive cross-section for production of any 
particle represented by Q in nucleus-nucleus collision and production of the same in nucleon-nucleon 
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collisions in the following form : 



where Q is the secondary, like pion/kaon/proton/antiproton, produced in high energy nucleus- 
nucleus(AB) or in proton-proton (PP) collisions. The term e{y,pT) could be expressed in the fac- 
torization form e{y,pT) = f{y)g{PT)- While investigating a specific nature of dependence of the 
two variables (y and py), either of them is assumed to remain averaged or with definite values. 
Speaking in clearer terms, if and when rapidity dependence is studied by experimental group, the 
transverse momentum is integrated over certain limits and is absorbed in the normalization factor. 
So the effective formula for rapidity spectra turns into 

-r\AB — >QX ~ {AB)f^y'^ — \pp — ^QX- (2) 
dy dy 

The main bulk of work, thus, converges to the making of an appropriate choice of form for f(y). 
And the necessary choice is to be made on the basis of certain premises and physical considerations 
which do not violate the canons of high energy interactions. 

Applying the concept of both limiting fragmentation and the Feynman Scaling hypothesis we 
attempt at providing first a fit to the data on rapidity distributions for proton-proton collisions at 
several y/s values, ranging from ^/s = 23 GeV to 63 GeV with a parametrization by 3-parameter 
formula : 

-— = Ci 1 + exp— — (3) 
a dy A 

where Ci is a normalization constant and yo, A are two parameters. The choice of the above form 
made by Thome' et al[4j was intended to describe conveniently the central plateau and the fall-off 
in the fragmentation region by means of the parameters yo and A respectively. For all five energies 
in PP collisions the value of A was obtained to be ~ 0.55 for pions[5j and kaons[6j, and ~ 0.35 for 
protons/antiprotons[6j. And these values of A are generally assumed to remain the same in the 
ISR ranges of energy. Still, for very high energies, and for direct fragmentation processes which 
are quite feasible in very high energy heavy nucleus-nucleus collisions, such parameter values do 
change somewhat prominently, though in most cases with marginal high energies, we have treated 
them as nearly constant. 

Now, the fits for the rapidity (pseudorapidity) spectra for non-pion secondaries produced in 
the PP reactions at various energies are phenomenologically obtained by De and Bhattacharyya[6] 



through the making of suitable choices of Ci and i/q. It is observed that for various secondaries the 
values of yo remain almost constant and do not show up any sharp species-dependence amongst 
the secondaries. However, for pions it gradually increases with energies and the energy-dependence 
of yo is empirically proposed to be expressed by the relationship |5] : 

yo = 0.55 In + 0.88 (4) 

However, the energy-dependence of yo is studied here just for gaining insights in their nature 
and for purposes of extrapolation to the various higher energies (in the frame of y^ijvTv) for several 
nucleon-nucleus and nucleus-nucleus collisions. The specific energy (in the cm. system, ^/snn) for 
every nucleon-nucleus or nucleus-nucleus collision is first worked out by converting the laboratory 
energy value(s) in the required cm. frame energy value(s). Thereafter the value of yo to be used 
for computations of inclusive cross-sections of nucleon-nucleon collisions at particular energies of 
interactions is extracted from Eq. (4) for corresponding obtained energies. This procedural step 
is followed for calculating the rapidity (pseudorapidity)-spectra for only the pions produced in 
nucleon-nucleus and nucleus-nucleus collisions [5]. But, for the studies on the rapidity-spectra of 
the non-pion secondaries produced in the same reactions one does neither have the opportunity to 
take recourse to such systematic step, nor could they actually resort to this rigorous procedure, as 
both the values of yo and A (as was done by De and Bhattacharyya[6J) do not depict any sensitive 
energy-dependences. So, we have considered only the average values of yo and A to study the 
rapidity spectra for non-pion secondaries produced in nucleon (nucleus)-nucleus collisions. 

Our next step is to explore the nature of /(y) which is envisaged to be given generally by a 
polynomial form noted below : 

/(y) = a + /3y + 7y2, (5) 

where a, /3 and 7 are the coefficients to be chosen separately for each AB collisions (and also for AA 
collisions when the projectile and the target are same). But, in so far as dealing with the rapidity 
spectra is concerned, the above expression suffers from a limitation which we would discuss later in 
some detail in section 3.2. Besides, some other points are to be made here. The suggested choice 
of form in expression (5) is not altogether fortuitous. In fact, we got the clue from one of the 
previous work by one of the authors (SB) [7] here pertaining to the studies on the behavior of the 
EMC effect related to the lepto- nuclear collisions. In the recent past Hwa et al[8j also made use of 
this sort of relation in a somewhat different context. Now we go back to our original discussion. 
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Combining Eqs. (2) and (5) the final working formula for ^ in various AB (or AA) collisions can 
be expressed by the following relation : 

^Ub^qx = C,{ABr^f^y^^y'^\pP^QX = CMBfy^^y\l + exp ^)-\ (6) 

where C2 is the normalization constant and C^=C2{AB)'^ is another constant as a is also a constant 
for a specific collision at a specific energy. The parameter values for different nucleus-nucleus 
collisions are given in Tablel - Table 9. However, in the next section (Section 3) we attempt at 
providing hints to the possible origin of our very starting equation [Eqn. 1] in this Section. 



3 Basic Sources of the Theoretical Framework : An OutUne 

At the very root, our approach owes a lot to the basics of the Glauber Model and to the pioneering 
developmental works by Gribov who showed how to incorporate the Glauber Model for interactions 
of hadrons with nuclei into the general framework of relativistic quantum theory. Subsequently, 
an important contribution to the theory multiparticle production was made by Gribov jointly with 
Abramovsky and Kancheli which is normally referred to as the AGK cutting rules. Later Capella 
et al|16) combined in an excellent work the outcomes of all of them in a concise paper from which 
we pick up the expressions having very close resemblance with both the physical ideas, the terms 
and also the mathematical structure. 

An important consequence of the space-time structure of the scattering diagrams for hadron- 
nucleus interactions as studied by Capella et al is : for inclusive cross sections all rescatterings 
cancel with each other and the cross sections are determined solely by the very common impulse 
approximation. This statement is valid asymptotically in the central rapidity region (y ~ 0) and 
for cases where the masses of the intermediate studies are both limited and energy-independent. 

The inclusive cross section for the production of a hadron i is expressed, for a given impact 
parameter b, in terms of inclusive cross section for hN interactions |12] 

where TA{h) is the nuclear profile function (J (PbTA{b) = A). After integrating over b we get 

cP'p d?p 
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The total and inelastic liA cross sections in the Glauber model can be easily calculated and are 
given for heavy nuclei by well known expressions. For example 

< = / dH{l - expi-aTMm (9) 

The situation for nucleus-nucleus collisions is much more complicated. There are no analytic 
expressions in the Glauber model for heavy-nuclei elastic scattering amplitudes. The problem 
stems from a complicated combinatorics and from the existence of dynamical correlations related to 
"loop diagrams" [13l m] . Thus, usually optical-type approximation [TSj 116] and probabilistic models 
for multiple rescatterings[17J are used. For inclusive cross sections in AB-collisions the result of 
the Glauber approximation is very simple to formulate due to the AGK cancelation theorem. It 
is possible to prove, for an arbitrary number of interactions of nucleons of both nuclei[18j, that 
all rescatterings cancel in the same way as for h A- interactions. Thus a natural generalization of 
eq. (7) for inclusive spectra of hadrons produced in the central rapidity region in nucleus-nucleus 
interactions takes place in the Glauber approximation 

where TAB{b) = J dP sTA{'^)TB{t — After integration over b eq. (10) reads 

E'i^ = ABE^^ (11) 

It is to be noted here that eqs. (10), (11) are valid for an arbitrary set of Glauber diagrams and 
are not influenced by the problem of summation of "loop" diagrams mentioned above. 

The densities of charged particles can be obtained from eqs. (10), (11) by dividing them by the 
total inelastic cross section of nucleus-nucleus interaction. For example 

dn%{h) ^ TAB{b)da'^^^ 

dy a% dy ^ ' 

and 

dn% ^ AB_da]^ 
dy a% dy ^ ' 

In the following we shall use these results to calculate particle densities in the central rapidity 
region at energies of RHIC and LHC. It is to be seen that our starting expression [eqn. (1)] is just 
a generalization of the expression (11) shown here. 
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4 Results 



This section is to be divided into three Subsections. In the first one we summarize the procedural 
steps that follow from the Section-2. And in the rest the results are delivered. 

4.1 

The procedural steps for arriving at the results could be summed up as follows : 

(i) We assume that the inclusive cross section (I.C.) of any particle in a nucleus-nucleus (AB) 
collision can be obtained from the production of the same in nucleon-nucleon collisions by multi- 
plying it (I.e.) a product of the atomic numbers of each of the colliding nuclei raised to a particular 
function, e{y,pT), which at first is unspecified (Equation 1). 

(ii) Secondly, we accept that factorization of the function e{y,pT) = f{y)g{pT) which helps us 
to perform the integral over px (Equation 2) in a relatively simpler manner. 

(iii) Thirdly, we assume a particular 3-parameter form for the pp cross section with the param- 
eters Ci, yo and A (Equation 3). 

(iv) Finally, we accept the ansatz that the function f(y) can be modeled by a quadratic function 
with the parameters a, /3 and 7 (Equation 5). 

4.2 

The results are shown here by the graphical plots with the accompanying tables for the parameter 
values. In Fig.l we draw the rapidity-density of the negative pions for symmetric Pb-|-Pb collisions 
at several energies which have been appropriately labeled at the top-right corner. In this context 
some comments are in order. Though the figure represents the case for production of negative 
pions, we do not anticipate and/or expect any strong charge-dependence of the results. Besides, 
the solid curves in all cases - almost without any exception - demonstrate our GCM-based results. 
Secondly, the data on rapidity-spectra for some high-energy collisions are, at times, available for 
both positive and negative y-values. This gives rise to a problem in our method. It is evident here 
in this work that we are concerned with only symmetric collisions wherein the colliding nuclei must 
be identical. But in our expression (6) the coefficient f3 multiplies a term which is proportional to 
y and so is not symmetric under y— >^(-y). In order to overcome this difficulty we would introduce 
here /3=0 for a large number of the graphical plots. These plots are represented by Fig.l, Fig.2, 
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Fig.5, Fig.6 and Fig.7 for the various secondaries produced in Au+Au/Pb+Pb collisions under 
different conditions of relatively low cm. energy values. The parameter values in this particular 
case are presented in the Tables 1,2,3,6,7,8,9. The Fig. 2a and 2b (both for /3=0) are for production 
of positive kaon and negative kaons in Pb+Pb interactions at 30A GeV and 20A GeV. But the Fig. 3 
is based on the data on pion-production in Au+Au collisions at several relatively low energies for 
which data are available on dn/dy (rapidity density of multiplicity) versus only positive y- values. 
The parameter values are presented in Table4. And the diagrams in Fig.4 represent the production 
of the major varieties of the secondaries produced in Au+Au collisions at RHIC; herein too the 
data are obtained for only positive values of rapidity. The corresponding parameter values are 
presented in Table5. 

4.3 

But the problem starts as soon as we attempt to reproduce the trends in data simultaneously for 
both +ve and -ve rapidity values. The graphs grow asymmetric even for symmetric collision. In 
Fig.5 with if='=-meson varieties we tried to fix and focus the problem encountered. The four labeled 
diagrams for ii'='=-meson varieties shown in the present Fig.5 provide fits to the data with unchanged 
A-values (A=0.55) for both positive and negative y-valucs quite separately. The parameter values 
are to be obtained from TablcG and TableT. These plots provide very important clues to our 
understanding the nature of A-values at very, very high energies. For both kaon-antikaons one 
observes quite clearly that the plots end up on the X-(rapidity)-axis for positive high-y- values 
within the range +5 to +6 (positive) y-values, whereas for the negative high y-values the spread 
widens to y-values ranging between -6 to -7. This shows up somewhat convincingly to us that 
the asymmetry problem arising in Fig.5 towards the -|-ve y-values cannot be addressed properly 
and remedied without inserting necessary changes in the accompanying sets of parameters, as the 
treated energy-ranges for these Au+Au collision cases are very high compared to the various sets 
of other data presented here. Our studied diagnosis of the root of the problem forced us to alter 
the previous A-values with phenomenological induction of increasing the nature of A with rising 
energy as was the case for yo- And, thus, making use of this ansatz, we were able to redeem the 
symmetry-feature of the rapidity plots. The subsequent two figures, the Fig.6 and Fig.7 strikingly 
illustrate how the qualitative improvement in the nature of the plots on rapidity-distributions for 
Au+Au reaction at RHIC was achieved by changing A-values. The parameter values for Fig.6 and 
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Fig.7 are given in TableS and TableQ. The Fig.6 and Fig.7 (for /3=0) describe the rapidity spectra 
on the various non-pion secondaries as labeled in Figs. 6(a), 6(b), 6(c), 6(d) and Figs. 7(a), 7(b), 
7(c), 7(d) for both positive and negative y- values in Au+Au collisions at RHIC energies with A- 
value=1.70 and A-value=3.5 respectively. Due to the paucity of data, we are not going to present, 
at this juncture, any quantitative prediction on the nature of rise of A, though we qualitatively 
accept and use it in this work. 

4.4 

The diagrams shown in Fig. 8(a) and Fig. 8(b) represent the model-based results on /■k'^ and 
K~ l'K~ respectively. These plots arc drawn on the basis of the figures shown in Fig. 4 with the 
fit-parameters given in Table5. The data-trends have been captured by them in a modestly right 
manner. The next plot (Fig.9) pertains to the particle/antiparticle ratio- values which are really 
quite crucial. This plot is drawn with the help of the figures shown in Fig.4 and the Table5, that 
is, they were worked out on the basis of the graphs drawn in Fig.4 and the accompanying table 
displayed in the Table5. 

The rest two plots. Fig. 10(a) and Fig. 10(b), are drawn on the basis of these revised figures with 
the gradually changed values of A. As no data on a charge-dependent basis for pions produced 
in Au-|-Au reactions at RHIC energies with negative y-values were available, we could not depict 
the 'K~ j-K^ ratio-behaviour in Fig.lO(a) and Fig.lO(b) with changed A-values. And the physical 
justifications for introduction of such enhancements in the values of A had already been hinted at 
quite emphatically in Section 2. They also reproduce the features of the data, in a fairly satisfactory 
manner. In fact, with changed A-values they describe, in our opinion, data-characteristics in a 
somewhat better way. 

5 Discussion and Conclusions 

Surveying the model-based plots and judging by the degree and nature of agreement between the 
measured data and the obtained results, we could fairly sum up by a reasonable statement that the 
model(s) applied here describe the data in a modestly satisfactory manner. The used models have 
both a deductive origin and some phenomenological components, with some conceptual physical 
strands : (i) the results for the nucleus-nucleus collisions could be arrived at by inducting the results 
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of nucleon-nucleon collisions and some relevant phenomenological product terms incorporating the 
physics of nuclear dependences; (ii) the property of factorization is also embedded in the results in 
an ingrained manner; (iii) there is a certain degree of universality of nature among the secondaries 
produced in nucleon-nucleon, nucleon- nucleus and nucleus-nucleus interactions at high energies, by 
virtue of which the entire ranges of rapidity-spectra have here been analyzed with the help of two 
parameters (yo and A) for PP collisions, and two other (/3 and y) parameters for nucleus-nucleus 
collisions. Of them, yo for pion-production cases alone could be obtained by a formula [egn.(4)] 
suggested by us. Similar behaviours of j/o and A are perfectly possible at very very high energies 
for all non-pion secondaries as well, though, in most cases, we leave them practically unchanged. 
However, There arc some exceptions for Fig. 6 and Fig. 7. Thus, we arc left with only two arbitrary 
parameters, which one should try to relate with the physics of nuclear geometry and collision 
dynamics. But there are some problems with the studies of the rapidity-spectra : firstly, the 
systematic studies on the rapidity spectra - both experimental and theoretical/phenomenological 
- are comparatively quite fewer in number ; secondly, the measurements suffer relatively from a 
higher degree of uncertainty. In fact, the centrality-dependent studies on rapidity spectra are not 
yet available, for which the studies on the relationship between nuclear geometry and the rapidity- 
spectra cannot yet be taken up either by us or others. Viewed from the angle of these existing 
limitations and constraints, we could humbly claim that the task of explaining the data-trends and 
behaviours on rapidity spectra is done here with a modest degree of success. We have checked that 
with the changed values of A, it is easier to describe the totality of data assembled here on several 
symmetric collisions at various energies within the domain of high energy ranges. 

In our work, we have categorically pointed out the problem of treating the cases of negative 
rapidity- values with the help of the generalized expression given by expression (6), wherein a 
specific term as is indicated clearly in 4.2 constitutes the limitation of the method. For some other 
similar cases, however, the method remains fully valid and works phenomenologically quite well. 

There is yet another drawback of this model. Since the coefficients are fit for different species 
and collision energies we cannot right now extrapolate to different collision energies or connect the 
various coefficients, nor could we make a reliable prediction for Pb-|-Pb collision at LHC on the 
basis of our fits to the RHIC data for Au-|-Au collisions. As a result, the model lacks, so far, any 
strong predictive power. And this is true of almost any phenomenological model whatsoever. But 
this difficulty could be remedied to a considerable extent, had there been availability of reliable and 
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accurate data at different higfi energies witfi reasonably spaced intervals for both PP interactions 
and/or any other symmetric high energy nuclear collisions. Besides, we are yet to find and work 
out the ways and means for adaptation of this methodology to any non-symmetric collision. 
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Table 1: Values of different parameters for production of negative-pions in Pb+Pb collisions at 
different energies(for (3 = 0) for both +ve and -ve rapidities. [Reference Fig. No.l] 



Energy {GeV) 


Cs 


7 


ndf 


6.3 


87.426 ±0.403 


-0.042 ±0.0006 


11.449/13 


7.6 


100.095 ±0.291 


-0.036 ±0.0003 


10.810/11 


8.7 


118.278 ±0.274 


-0.036 ±0.0002 


2.943/5 


12.3 


147.107± 0.401 


-0.026 ±0.0002 


10.218/6 


17.3 


184.850 ±0.449 


-0.023 ±0.0001 


3.637/4 



Tabic 2: Values of different parameters for production of kaons in Pb+Pb collisions at different 
energies for both +ve and -ve rapidities (for /?=0) . [Reference Fig. No. 2(a)] 



Energy 


Cz 


7 


xL 

ndf 


2QAGeV 


5.798 ±0.015 


-0.0903 ±0.0006 


0.673/3 


^QAGeV 


8.000 ±0.027 


-0.0726 ±0.0005 


4.225/5 



Table 3: Values of different parameters for production of anti-kaons in Pb+Pb collisions at different 
energies for both +ve and -ve rapidities (for /?=0). [Reference Fig. No. 2(b)] 



Energy 


Gs 


7 


•2 

ndf 


20AGeV 


16.990+0.059 


-0.0430 ± 0.0020 


6.888/8 


SOAGeV 


21.589 ±0.066 


-0.0404 ±0.0015 


7.282/11 



Tabic 4: Values of different parameters for production of negative-pions in Au+Au collisions at 
different energies for +ve rapidities alone (for /? 7^ 0). [Reference Fig. No. 3] 



Energy{GeV) 


C3 


13 


7 


xi 

ndf 


8.8 


110.243+0.184 


0.013 + 0.0002 


-0.032 + 0.0002 


2.676/11 


12.4 


147.023+ 0.369 


0.005 + 0.0002 


-0.019 + 0.0001 


6.494/9 


17.3 


185.589+ 0.252 


-0.006+0.0001 


-0.009 + 0.0008 


3.238/8 



Tabic 5: Values of the different parameters for production of pions, kaons, proton-antiprotons in 
Au+Au collisions at ^/Snn = 200GeV for +ve rapidities alone (for /3/0). [Reference Fig. No.4] 



Production 


Cz 


/3 


7 


xi 

ndf 


7r+ 


303.480+0.400 


0.0009 + 0.00007 


0.0052 + 0.00002 


0.301/4 


7r~ 


300.938+ 1.981 


0.0009 + 0.0003 


0.0057 + 0.0001 


4.668/4 


K+ 


49.169+0.099 


0.0011 + 0.00009 


0.0049 + 0.00003 


0.073/4 


K- 


46.996 + 0.073 


0.0010 + 0.00009 


0.0027+0.00003 


0.026/3 


P 


24.862 + 0.659 


0.009 + 0.002 


0.0059 + 0.001 


3.431/8 


P 


17.794 + 0.290 


0.011 + 0.003 


-0.0059 + 0.001 


0.667/6 
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Table 6: Values of different parameters for production of kaon and antikaon in Au+Au collisions 
at -s/Sjvjv = 200GeV (for j3 = 0) for +ve rapidity only, taking A=0. 55. [Reference Fig. No.5] 



Production 


Co. 


7 


ndf 


Kaon 


47.128 iO.Ool 


-0.005 ± 0.0001 


1.941/09 


Anti — kaon 


45.207 ±0.064 


-0.008 ±0.0001 


1.018/06 



Table 7: Values of different parameters for production of kaon and antikaon in Au+Au collisions 
at \/Snn = 200GeV (for (5 = 0) for -ve rapidity only, taking A=0. 55. [Reference Fig. No.5] 



Production 


Cs 


7 


xi 

ndf 


Kaon 


47.808 ±0.076 


-0.008 ±0.0001 


1.990/09 


Anti — kaon 


46.705 ±0.069 


-0.011 ±0.0001 


3.461/08 



Table 8: Values of different parameters for production of proton, antiproton, kaon and antikaon 
in Au+Au collisions at \/Snn = 200GeV (for ^ = 0) for both +ve and -ve rapidities, taking 
A= 1.7. [Reference Fig. No.6] 



Production 


Cs 


7 


xi 

ndf 


Proton 


49.770 + 0.690 


-0.011 ±0.002 


2.083/09 


Anti — proton 


21.103 ±0.408 


-0.016 ±0.001 


4.636/13 


Kaon 


55. 524± 0.111 


-0.009 ±0.0003 


2.867/04 


Anti — kaon 


48.660 ±0.109 


-0.010 ±0.0003 


3.859/05 



Table 9: Values of different parameters for production of proton, antiproton, kaon and antikaon 
in Au+Au collisions at ^/Snn = 200GeV (for ^ = 0) for both +ve and -ve rapidities, taking 
A=3. 5. [Reference Fig. No. 7] 



Production 


Cs 


7 


ndf 


Proton 


59.529+ 1.028 


-0.009 + 0.001 


4.924/13 


Anti — proton 


25.107 + 0.469 


-0.016 + 0.001 


5.160/13 


Kaon 


67.962 + 0.228 


-0.011 + 0.0003 


18.708/09 


Anti — kaon 


66.694 + 0.198 


-0.013 + 0.0003 


9.524/05 
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Figure 1: Rapidity spectra for secondary pions produced in Pb+Pb interactions at different energies 
for j3 =0. The different experimental points are from Ref.[19j and the parameter values are taken 
from Tablel. The solid curves provide the GCM-based results. 



Pb+Pb Collisions 
30AGeV*1.4 

20AGeV ^-»— 
GCM 





(a) 



(b) 



Figure 2: Rapidity distributions of and produced in central Pb+Pb collisions at 20A and 
30A GeV(both for /? = 0). The different experimental points are from Ref. [20j and the parameter 
values are taken from Table2 and TableS. The solid curves provide the GCM-based results. 
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Figure 3: Rapidity spectra for negative secondary pions produced in Au+Au interactions at different 
energies (for The different experimental points are from Ref.|2.1j and the parameter values 

are taken from Table4. The solid curves show the GCM-based results. 
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Figure 4: Plot of vs. y for secondary pfon, kaon, proton and antiproton produced in Au+Au 
interactions at \/Snn = 200GeV (for /3y^0). The different experimental points are from Ref . |22j 
and the parameter values are taken from Table5. The solid curves provide the GCM-based results. 
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Figure 5: Plot of vs. y for secondary kaon and antikaon produced in Au+Au interactions for 
positive and negative y- values separately at \/Snn = 200GeV (for (3=0), taking A=0.55. The 
different experimental points are from Ref. |23j and the parameter values are taken from Table6 
and Table?. The solid curves provide the GCM-based results. 
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Figure 6: Plot of ^ vs. y for secondary proton, antiproton, kaon and antikaon produced in Au+Au 
interactions for positive and negative y- values separately at \/Snn = 200GeV (for /3=0), taking 
A=1.7. The different experimental points are from Ref. \23\ [M] and the parameter values are taken 
from Tables. The solid curves provide the GCM-based results. 
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Figure 7: Plot of vs. y for secondary proton, antiproton, kaon and antikaon produced in Au+Au 
interactions for positive and negative y- values separately at \/Snn = 200GeV (for /3=0), taking 
A=3.5. The different experimental points are from Ref. \23\ [M] and the parameter values are taken 
from Table9. The solid curves provide the GCM-based results. 
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Figure 8: Full phase-space K/tt ratios as a function of rapidity systematics at ^JSnn = 200GeV. 
Errors are statistical only. The different experimental points are from Ref.[25] and the parameter 
values are taken from Table5. The solid curves provide the GCM-based results. 





Figure 9: Antiparticle-to-particle ratios as a function of rapidity in central Au+Au collisions at 
\/Sj\fN = 200 GeV. The solid lines are our results for 7r^/7r+, /K^, Pbar/P respectively. The 
experimental data are from Ref . |22) and the parameter values are taken from Table5. 
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Figure 10: Antiparticle-to-particle ratios as a function of rapidity in central Au+Au collisions at 
VSnn = 200 GeV. The solid lines are our results for K~ /K'^ and Pbar/P. The experimental data 
are from Ref. [22| and the parameter values are taken from TableS and TableO, taking A=1.7 and 
A=3.5 for Fig. 10(a) and Fig. 10(b) respectively. 
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